s 



PUPT-1479 
hep-th/94070T4 



A modified c = 1 matrix model 
^ : with new critical behavior 



^^ '. Steven S. Gubser and Igor R. Klebanov 

^ ! Joseph Henry Laboratories 

"^ ! Princeton University 

O ' Princeton, New Jersey 08544 

O 



^ ■ Abstract 



k> ! By introducing a /(it (7 (Tr $^(t)) term into the action of the c = 1 matrix model of 

two-dimensional quantum gravity, we find a new critical behavior for random surfaces. The 
planar limit of the path integral generates multiple spherical "bubbles" which touch one 
another at single points. At a special value of g, the sum over connected surfaces behaves 
as A^log A, where A is the cosmological constant (the sum over surfaces of area A goes as 
A~^). For comparison, in the conventional c = 1 model the sum over planar surfaces behaves 
as AVlogA. 
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1 Introduction 

In recent years a considerable effort has been devoted to understanding the statistical me- 
chanics of random surfaces. This problem is relevant to Polyakov string theory in non-critical 
dimensions. Thanks to the matrix model techniques, we now have a thorough understanding 
of random surfaces embedded in one dimension [|I], |^, |^ or less 0]. Little is known, however, 
about the physically more interesting higher dimensional embeddings. With this in mind, it 
is important to continue formulating and solving new matrix models. 

One interesting modification of the conventional one- matrix model was solved in ref. [^. 
A new type of critical behavior arises when a term of the form (7(Tr$^)^ is added to the 
action, leading to an integral over the N x N hermitian matrix $ of the form 



-N 



Tr(i#2-iA$4)_|,(Tr<I.2) 



log /P$e I y^ '^ ^ ^y n . (1) 

The large A^ limit of this quantity has an interesting geometrical interpretation. Feynman 
graphs of the perturbation theory in A generate the usual connected closed random surfaces, 
while the g{TT^^y term can glue a pair of such surfaces together at a point. This point 
can be resolved into a tiny neck (a wormhole), so that the network of such touching surfaces 
can be assigned an overall genus. In the leading large A^ limit one picks out the surfaces of 
genus zero, which look like trees of spherical bubbles such that any two bubbles touch at 
most once, and a bubble is not allowed to touch itself. 

The authors of ref. 0] found a critical line in the (A, g) plane where the free energy Eq. (|l|) 
becomes singular. There exists a critical value gt such that, for g < gt, the singularity in the 
function of A is characterized by 7str = — 1/2. In this phase the touching of random surfaces 
is irrelevant and one finds the conventional c = behavior. For g > gt, on the other hand, 
7str = 1/2, and one finds branched polymer behavior, which is dominated by the touching. 
Most interestingly, for g = gt, the authors of ref. P found a new type of critical behavior 
with 7str = 1/3. The interpretation of this is not completely clear, but it has been suggested 
that at this point one has an effective theory for random surfaces embedded in more than 
one dimension ^. In fact, after some fine tuning, theories with 7str = 1/"^ can be formulated 

In view of these results, a continued study of the f7(Tr$^)^ term is warranted. In this 
paper we investigate its effect on random surfaces embedded in one dimension. The relevant 
model is the matrix quantum mechanics defined by the path integral 



where $ is an A^ x A^ hermitian matrix. T will be taken to oo in the end, but we retain 
it as a finite quantity for the time being. In the large A^ limit we find a critical line where 
the free energy is singular. As in the model Eq. (|l|), there exists a critical value gt which 
separates the conventional matrix model behavior (in this case c = 1) from the branched 



polymer phase. For g = gt we find a new critical behavior: the sum over connected surfaces 
of genus zero behaves as A^logA, where A is the cosmological constant. This should be 
compared with the usual c = 1 behavior, A^/logA. Our modified critical behavior is much 
simpler after we carry out the inverse Laplace transform: the sum over surfaces of fixed 
area A behaves as A~^. The corresponding formula for c = 1 is more complicated due to 
logarithmic corrections. 

The organization of this paper is as follows. In section ^, we show how the modified c = 1 
model can be reduced to a fermionic system. In section ^ we find conditions on the ground 
state of the fermionic hamiltonian which allow us to find the critical points of the model and 
evaluate 7str- We end with a brief discussion in section §. 

2 Reduction to fermions 

The sum over connected surfaces is given, up to a factor, by the free energy 

y'=^\ogZ. (3) 

In the N ^ oo limit, the leading term of JF (which is the only term we will consider in this 
paper) is 0{N'^). This dominant term generates touching spherical bubbles embedded in one 
dimension. To be more explicit, the single trace terms in Eq. (^) generate planar, quartic 
Feynman graphs embedded in one dimension: these are the individual bubbles, which look 
just like ordinary c = 1 surfaces. The double trace term, (7(Tr$^) , gives rise to a pair 
of "touching" propagators. When these propagators are incorporated into two different 
bubbles, we interpret the bubbles as touching. Arrangements where a pair of bubbles touch 
in more than one place, or where a bubble reaches around and touches itself, are suppressed 
by a factor l/N"^ and thus correspond to 0(1) corrections to JF. 

The path integral in Eq. (H) can be written as a transition amplitude 

Z = (/|e-^^^|^) where 



H = L-^ + Tr (l^^-lx^'\-^(i:i^^y . (4) 

In the T —>■ oo limit, JF is — A^ times the ground state energy of this hamiltonian. The ground 
state must be a S'[/(A^)-invariant function of $, which is to say a symmetric function of the 
eigenvalues Xj of $: dependence on "angular" degrees of freedom can only raise the energy. 
As in the c = 1 model, the key step is to pass to a fermionic system, in which the ground 
state is an anti-symmetric function of the Xj. If A(xi) is the Vandermonde determinant of 
the eigenvalues Xi of $, then 



where Hf is the fermionic hainiltonian: 

The (Tr$^) term introduces interactions among the fermions. These interactions can 
be taken care of by a self-consistent field approach, analogous to Hartree-Fock calculations 
in multi-electron atoms. For the purpose of finding the ground state energy of Hf to leading 
order in A^, it is permissible to make the replacement 

(Tr$2)^^2(Tr<l>2)Tr$^- (Tr<l>2)2 . (7) 

An intuitive way to justify this replacement is to consider Tr $^ as varying slightly around 
its expectation value: Tr$^ = (Tr$^) +(5Tr$^. Expanding (Tr$^) to first order in (5Tr $^ 
yields Eq. (0). 

Applying Eq. (J^) to Eq. @ turns Hf into a constant plus a sum of A^ single-particle 
hamiltonians: 
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Hf-^gN( — Tr $2 ) + V -— ^ + -x'j - -Xxf - 2g ( — Tr $ 
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We write the single-particle hamiltonian as 

Mx) = -^^ + ^(1-4^c)x2-1ax^ where c=(lTr$2). (9) 

Since our modified matrix model reduces to free fermions, its solution proceeds essentially 
along the lines of the usual c = 1 solution. However, the necessity of imposing a self- 
consistency condition adds an extra ingredient which, as we will show, can modify the nature 
of the critical behavior. 

3 Leading order solution 

To obtain the ground state energy of Hf to leading order in A^ as A^ ^ 00, it suffices to use 
semi-classical techniques. We regard the A^ fermions with single-particle hamiltonian h{x) 
as a Fermi gas in the potential U{x) = |(1 — Agc)x^ — \\x'^. If the Fermi energy is ep, then 
the particle density in x is 

A^ I 

pW = -V2(eF - U{x)) . (10) 

Given A and g, the values of ep and c are determined by a normalization condition and a 
self-consistency condition: 



1 /■ 1 

— / dx p{x) = 1 and — / dxx'^p{x) = c . (11) 



If we let E be 1/A^ times the ground state energy of Hf, then 

-L. = E = gc' + eF-^jdx [2(e^ - U{x))f' . 
As a check on the self-consistent field method, one may easily show that 



(12) 
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— / dxx^p(x) = — ( — Tr $^ 

NJ '^^ ' 4:\N 
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(13) 



These results agree with what one finds by explicit differentiation of Eq. (0) with respect to 
A and g. 

At this point it is convenient to rescale x and h{x) as follows: 
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+ V{y) where 



(14) 



A 



and V{y) = -y^ - -y^ . 



We also define a rescaled Fermi energy jip = (i_4 c)^ ^f- The highest possible value for fxp is 
fie = 1/4: at this energy, the fermions completely fill the local well of V{y). 
Upon rescaling, Eq. (|ll|) becomes 



A 



1 rv' 



(1 - %c)3/2 
A^c 
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dyd2{fip-V{y)) = I{fiF) 



IX J-y, 



dyyy2{fiF-V{y)) = J{fiF) 



where yc 



1 — y/1 — A/jp. c can be eliminated from these equations to yield 
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(15) 



(16) 



Eqs. (|15D and (piq) represent necessary conditions on the ground state, but as we shall see, 
they are not quite sufficient, as there are sometimes two possible values of fip for a given A 
and g. It will not be hard to tell which represents the true ground state. 

Eq. (p!6D determines a family of curves in \-g space that depends on the single parameter 
fip- Fig. d^) shows twenty of these curves. The outer envelope of all the curves is the locus 



of critical points where the perturbation series that generates the random surfaces begins to 
diverge. To see this, suppose one varies A and g from toward criticahty. As long as the 
point (A, (?) is inside the outer envelope, there is a ground state of if/, and E is finite. But 
as soon as (A, g) crosses the envelope, there is no solution to Eq. ([T5|), so E and hence JF are 
undefined. 

There are two regions of the critical curve. For g < gt = j^ ( |j , the curve given by 
Eq. (0) with fip = fic is the critical curve. For g > gt, the critical curve is tangent to a 
curve given by Eq. ([iBD with ftp < fic- Because of the different critical behaviors we shall 
find in these two regions, we shall call the region g < gt the c = 1 region and the region 
g > gt the branched polymer region. While the critical behaviors in these two regions have 
been found in other models, right at g = gt we find something new. 

In all the calculations we present, we continue to use the quartic potential, ATr<l>^, that 
we started out with in Eq. (H). As a check on universality, however, we have performed the 
same calculations with a cubic potential, and we found the same regions of the critical curve 
and the same three types of critical behavior. 

3.1 Critical behavior in the c = 1 region 

Let us pick a point (Acf^c) on the c = 1 part of the critical curve. To determine the critical 
behavior of £■ as A — ;► Ac with g fixed at gc, it suffices to calculate dE/dX, which from 
Eq. (|13D can be shown to be 



— = --y - where K = ir(/i^) ^-j_Jyy'p{^p - V{y)) . (17) 
To analyze the behavior of /, J, and K for ^p near /ic, let us set fip = f^c — fJ' and write 



Ic = /(^e) = ¥ , Jc= J(/i.) = '£, K, = Kif^c) = i? , ^^g^ 

I{fic - n) = Ic + SI , J{fic - ^) = Jc + SJ , K{iic - fi) = Kc + 6K . 
A standard result from c = 1 calculations is 

dl 1 fV'^ dy 1 , ^, , , , 

7^ = / , = = ^log/x + 1 , 19 

which means 61 = -^r/^filogfi + 0{ji). To calculate 5J and 5K, let us first define yc = 
1 + VI - 4/iF, so that fiF - V{y) = \{y^, - y^)ig - y^). Now, 

dJ dl V2 p {l-y^)dy 



-^ — / dy = 3/c as /i ^ 
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- ^C dy{l + y')=Ah 
Keeping terms of order /x and larger, we find 
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(20) 



5J = 5I + 3/c/i 



and 



5K = 5I + 4J,/i 



(21) 



Now we perform a perturbation of A away from Ac with g fixed at g^- If one defines 
z = (A//)^/3 and Zc = {Xc/Icf'^, Eq. (|l|) becomes 



9j2 =z{l-z^) 



(22) 



Using this along with Eq. ([21| ) one can derive 

5z 



- = ^~^' (^3^ + 15 ^ 

Zc 1 - 3^2 I /^ , 



5A ^5 -6^2 5/ ,r^-zl 

V = 2^ :ri-r + 45- -^^i 

Xr 1 — 32;;? /c 1 — 6zi 



(23) 



In this and all the rest of our calculations in this section, we retain terms only up to 0{fi). 
Before continuing with the calculation of the critical behavior of E, let us note how 
Eq. (p3D can be used to determine the point {Xt,gt) where the fip = /^c curve ceases to 
coincide with the critical line. Eq. (|23|) actually applies to any point (Ac, gc) on the fip = fie 
curve, except when 1 — Sz"^ = 0. If we set 6X = 0, then Eq. (p3D becomes 







2- —^— + 45: 



/^ 



(24) 
l-32c"/c ■ l-32c'' ' 

which is a condition on Ac that determines where the fip = fJ^c — fJ^ curve intersects the 

fip = fic curve. Now, (Aj, gt) is the point where the fip = fie curve just starts to meet curves 

with lower fip. Hence A^ can be determined by letting /i — >• in Eq. (^41). Since the dominant 

term on the right hand side of Eq. (|^) is the one involving 61 /Ic, we must have z"^ -^ 5/6, 

whence A* = 3^ (|) and gt = jt (f) • 

Now let us see how E behaves near criticality. Retaining terms up to 0{fi) as usual. 
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Eqs. ( PB| ) and ( PSJ ) lead to the existence of two different types of critical behavior. First, let 
us take Ac > Af {gc < gt)- In this case Eqs. (EBI) and (U^) combine to give 
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dX AzJ. /3 \ 3 h-Qzf. X^ 3 5-6^2' 

As in the case of ordinary c = 1, we define the cosmological constant A by A = Ac — A. To 
leading order in A, 

_ n 1 1 - 3z^^ A 
^~~72if5-6z2 1ogA ' ^^'^^ 

which is the behavior typical of ordinary c = 1. The leading nonanalytic behavior of E is 

^ / 1 .^ QtT^ 1 fl-3zjy A2 ^, . , 

i? = (analytic m A) ^— r r- + (higher order terms) , (28) 

32\/2 z^ \5 - 6z^ J log A 

which again is typical of c = 1. Our conclusion is that, for g < gt, the touching of random 
surfaces is irrelevant in the sense that it does not destroy the c = 1 behavior. 

If we fix (yf = gt and let A — > A^, then we get a different critical behavior: Eq. (^3|) becomes 
6X/Xt = —5fj,, so that the leading order relation between A and fj, is just 

The fact that A and fi are simply proportional to each other is the crucial feature of the 
new critical behavior. In contrast to Eq. (^), the universal part of E is now 

36 /3\V2 
E = (analytic in A) H — g-^ ( - J tt^ A^ log A + (higher order terms) . (30) 

Similar calculations are possible for perturbations of g away from criticality with A held 
fixed. If one sets g = gc — T and A = Ac, then in the c = 1 region one finds that the singularity 
of E is rVlogr. For A = Ai, E ~ r^ logT instead. 

3.2 The branched polymer region 

It is clear from Fig. (|T[) that the region above the curve given by Eq. (|T^) with fj,p = fie, but 
below the critical curve, is covered not once but twice by the family of curves in Eq. (p!6|). 
Given (A, g) in this region, there are thus two values oi fip and c that satisfy the self-consistent 
field equations Eq. (|T^. The true ground state of H corresponds to the solution with the 
lower energy E. It is an unsurprising but slightly non-trivial fact that the solution with 
lower E is the one with lower fip. At any given point in Fig. (|^) where two curves cross 
each other, the more down-sloping curve is the one with lower np, and thus represents the 
ground state. 

Let us fix g and consider how A varies with fip. The singularity of the perturbation 
theory in A occurs at the maximum of X{fip). For any g, A(0) = 0. The behavior of A with 
increasing fip, however, depends on the value of g. 



For g < Qt, X is an increasing function of fip and the maximum occurs at the end-point 
fJ'F = fJ'c- The slope diverges logarithmically as fip approaches /ic, which leads to the c = 1 
critical behavior. For g = gt, A is an increasing function of fip with a slope that is finite 
everywhere. The finiteness of the slope at fip = fJ-c leads to the modified critical behavior. 

For g > gt, X^fip) increases to a quadratic maximum at a value of fip less than /Xc- This 
maximum corresponds to A reaching the critical curve: X{fim) = Xc- If we increase hf past 
fim, A(/ii?) decreases; this region represents only unphysical solutions to the self-consistent 
field equations. In other words, if for a given value of A there are two possible values of ftp, 
then we pick the smaller one because it is analytically connected to the perturbative region 
near A = 0. A graph of X^fip) with g = ^gt is shown in Fig. (|). 

The main point is that we can perturb A away from its critical value Ac by lowering fip: 
if we write A = Ac — A and fip = fJ^m — /^, then 

A = a/i^ -|- (higher order terms) , (31) 

where a > is some constant. 

From Eq. ( PT| ) and Eq. (^Tj) it is simple to show that 



dE IK 
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^Ja dfip P/^ 



A + 0(A) . (32) 



K/P^^ always has positive derivative, so the universal behavior of E is 

E = (analytic in A) — aA^'^ + (higher order terms) , (33) 

where a is a positive number. Predictably, the same power law occurs in the universal part 
of E when we perturb g = gc — T with A fixed at Ac: the leading nonanalytic behavior of E 

is ~ -r3/2. 

It is remarkable that as soon as gc exceeds gt, the Fermi level fip starts falling. Since the 
universal behavior of the c = 1 model is driven by the Fermi energy approaching criticality, 
we regard the dropping of the Fermi energy as the reason for the transition to the branched 
polymer phase. 

4 Discussion 

By performing the inverse Laplace transform of Eq. (^) we find that, for g = gt, the sum 
over connected genus zero surfaces of area A is 

This is precisely the KPZ scaling p[ with 7str = 0. For the conventional c = 1 matrix model 
we instead have M 



•^(^) ~ ,,,/ ,^o • (35) 

Thus, the critical behavior we have found is actually simpler than the conventional c = 1 
behavior. Is it possible that we have penetrated the c = 1 barrier, as suggested in ref. 0? 
Although we do not have a definite answer, it seems likely to us that the new critical point 
can be described by two-dimensional string theory. The scaling violations for c = 1 have 
been attributed to the unusual dependence of the tachyon potential on the Liouville 
field, T((/)) ~ (\)e'^ . It appears that, if the tachyon potential has the ordinary Liouville form. 



T(0) ~ e''^, then the simpler scaling of Eq. (^J) follows. Clearly, more work is needed to 



settle the stringy interpretation of our new critical theory. 

We believe that there are more interesting calculations that can be performed in the 
modified c = 1 matrix model. Such calculations should shed light on the critical behavior 
for higher genus surfaces and, ultimately, on the double-scaling limit. 

Note added 

After this paper was completed, we learned of an interesting paper by F. Sugino and 
O. Tsuchiya |10| in which results similar to ours were obtained by means of collective field 



theory. 
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Figures 

1. The family of curves representing all possible solutions {X,g) to the self-consistent 
field equations. Each curve has constant Fermi energy hf- The critical curve is the 
outer envelope of all these curves. The transition point from the c = 1 region to 
the branched polymer region is {\t, 9t) = ( 3^ (|) , j^ (|) )) where curves with 
different fip start crossing one another. When two curves do cross, the one with the 
lower fip (or, equivalently, more negative slope) represents the physical ground state. 
Thus, for each curve, only the part to the right of the contact with the envelope is 
physical. 

2. A graph of A versus fip for g fixed at ^gt- Only the points with fj,p < Hm represent 
physical solutions to the self-consistent field equations. 
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